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Masses and decay constants of the scalar quarkonia, XQoiQ ~ c) with quantum numbers 
I'^{J^'~") = 0^(0"''^) are calculated in the framework of the QCD sum rules approach both in vacuum 
and finite temperature. The masses and decay constants remain unchanged up to T ~ 100 MeV 
but they start to diminish with increasing the temperature after this point. At near the critic or 
deconfinement temperature, the decay constants reach approximately to 25% of their values in vac- 
uum, while the masses are decreased about 6% and 23% for bottom and charm cases, respectively. 
The results at zero temperature are in a good consistency with the existing experimental values 
and predictions of the other nonperturbative approaches. Our predictions on the decay constants in 
vacuum as well as the behavior of the masses and decay constants with respect to the temperature 
can be checked in the future experiments. 

PACS numbers: 11.55.Hx, 14.40.Pq, ll.lO.Wx 
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I. INTRODUCTION 

To understand the internal structure of scalar mesons has been a prominent topic in the last 30-40 years. Although 
the scalar mesons have been investigated for several decades, many properties of them are not so clear yet and 
identifying the scalar mesons is difficult, experimentally. Hence, the theoretical works can play a crucial role in this 
respect. In particle physics, the quarkonia refers to flavorless mesons containing a heavy b (c) quark and its own 
antiquark, i.e., bb (bottomonium) and cc (charmonium) . These approximately non-relativistic systems are the best 
candidates to investigate the hadronic dynamics and study the perturbative and non-perturbative aspects of QCD. 
It was believed that the quarkonia can help us to extract the nature of quark-antiquark interaction at the hadronic 
scale and play the same role in probing the QCD as the hydrogen atom play in the atomic physics 

A large number of the beauty and charmed systems have been experimentally observed in the last few decades (see 
for instance the theoretical calculations on the properties of these systems have been made mainly using 

potential model, where the quarkonia is described by a static potential, V = + kr and its extensions like the 

Coulomb gauge model 1349] • The first term in the potential is related to one gluon exchange and the second term is 
called the confinement potential. The recent CLEO measurements on the two-photon decay rates of the even-parity, 
P-wave scalar 0+"*", Xb{c)o ^-nd tensor 2++, Xb(c)2 states ([IS El and references therein) were motivation to investigate 
the properties of the quarkonia and their radiative decays from the quark-antiquark interaction point of view (see for 
example 

In [l^, which is a recent study on extraction of ground-state decay constant from both sum rules and potential 
models, it is stated that results obtained at each step of the extraction procedure both in QCD and in potential models 
follow the same pattern, hence all our findings concerning the extraction of bound-state parameters from correlation 
functions obtained in potential model can apply also to QCD. It is also proven that in QCD sum rules approach 
by tunning the continuum threshold which is related to the energy of the first exited state specially with a Borel- 
parameter dependent threshold, we can get a more reliable and accurate determination of bound-state characteristics 
comparing the potential models. The QCD sum rules approach as a non-perturbative approach is one of the most 
powerful and applicable tools to spectroscopy of hadrons and can play a crucial role in investi gati on of the properties 
of the hadrons jlSl - llq . It has been used to calculate the masses and decay constants of mesons jl9l - t26| . This approach 
was extended to contain the properties of the hadrons at finite temperature called thermal QCD sum rules |27H29l | 
supposing that the operator product expansion (OPE) and the quark-hadron duality assumption remain valid, however 
the quark-quark, quark-gluon and gluon-gluon condensates are altered by their thermal versions. The main aspiration 
of this addendum was to explain the results obtained from the heavy ion collision experiments. It is presently believed 
that the hot and dense medium where the hadrons are formed modifies masses and decay widths of hadrons. It 
is shown that heavy mesons like J/tp and also radial and orbital cc excitations have different behavior when the 
temperature of the medium changes (see [s^ and references therein). In [3l|, scalar mesons and scalar glueballs are 
investigated in holographic QCD at finite temperature . A flood of papers have also been dedicated generally to the 
determination of the condensates, mass and decay constant of mesons and some properties of the nucleons at finite 
temperature (H-H^. 

In the present work, we calculate the mass and decay constant of the heavy scalar XQO mesons with quantum 
numbers l'^{J^'~^) = 0"'"(0"'""'") using the thermal QCD sum rules approach. Here, we assume that with replacing the 
vacuum condensates and also the continuum threshold by their thermal version, the sum rules for the observables 
(masses and decay constants ) remain valid. In calculations, we take into account the additional operators in the 
Wilson expansion at finite temperature |46j and modify spectral density in QCD side. These operators are due to 
the breakdown of Lorentz invariance at finite temperature by the selection of the thermal rest frame, where matter is 
at rest at a definite temperature (ssl. l47j. In this condition, the residual 0(3) invariance brings these extra operators 
with the same mass dimension as the vacuum condensates. We also consider the interaction of the currents with the 
existing particles in the medium at finite temperature. Such interactions require modification of the hadron spectral 
density. 

The outline of the paper is as follows: in next section, sum rules for the the mass and the decay constant of the 
heavy scalar, XQo niesons are obtained in the framework of the QCD sum rules at finite temperature. Section HI 
encompasses our numerical predictions for the mass and decay constants as well as comparison of the results with the 
existing predictions of the other non-perturbative approaches and experimental values. 



II. QCD SUM RULES FOR THE MASS AND DECAY CONSTANT 

In this section, we obtain sum rules for the mass as well as the decay constant of the scalar quarkonia containing 
b or c quark in the framework of the thermal QCD sum rules. For this aim, we will evaluate the two-point thermal 
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correlation function, 

n(g,T) = z / d4xe'«-(r(J^(x)J^(0))), (1) 



in two different ways: pliysical and theoretical representations. In correlation function T denotes the temperature, 
T is the time ordering product and J^{x) = Q{x)Q(x) is the interpolating current of the heavy scalar meson, 
S — XQoiQ = ^1 c). The thermal average of any operator, A can be expressed as: 

where H is the QCD Hamiltonian, and (3 — 1/T is the inverse of the temperature T and traces are carried out over 
any complete set of states. 

The physical or phenomenological representation of the aforementioned two-point correlation function is obtained 
in terms of the hadronic parameters saturating it with a tower of scalar mesons with the same quantum numbers as 
the interpolating current. The theoretical or QCD representation is gained via operator product expansion (OPE) in 
terms of the QCD parameters such as quark's masses, and the vacuum condensates considering the internal structure 
of these mesons, i.e., quarks, gluons and their interactions with each other as well as with the QCD vacuum. Sum 
rules for the physical observables such as the decay constant and mass are obtained equating these two different 
representations through dispersion relation. To suppress the contribution of the higher states and continuum, Borel 
transformation with respect to the Qq = —q^ is applied to both sides of the sum rules for physical quantities. 

To calculate the phenomenological part, we insert a complete set of intermediate states owing the same quantum 
numbers with current between the currents in Eq. ([T]) and perform the integral over. As a result, at T = 0, we 
obtain 

n(,,o)^<»'^WI^)(^|J(o)|o)^...^ (3) 

vrig — q 

where • • • represents the contributions of the higher states and continuum and nis is mass of the heavy scalar meson. 
The matrix element creating the scalar meson from the vacuum can be written in terms of the decay constant, fs by 
the following manner: 

(0 I J(0) I S) = fsms. (4) 

Note that Eqs. ([3]) and Q are valid also at finite temperature, hence, the final representation for the physical side 
can be written in terms of the temperature dependent mass and decay constant as: 

n(,r, = Mm,... ,5) 

In QCD side, the correlation function is calculated in deep Euclidean region, q^ <C —^qcd '^i^ OPE where the 
short or perturbative and long distance or non-perturbative effects are separated, i.e., 

uQ'^D{q,T) = nP'=''*(g,T) +n"°"P'=''*(g,T). (6) 

The short distance contribution (bare loop diagram in figure ([1]) part (a)) is calculated using the perturbation theory, 
whereas the long distance contributions (diagrams shown in figure ([1]) part (b) ) are represented in terms of the 
thermal expectation values of some operators. To proceed, we write the perturbative part in terms of a dispersion 
integral, 

uQCD^q^ T) = / ^ffifill + n"°"P-'-*, (7) 
J s- q^ 

where, p{s,T) is called the spectral density at finite temperature. The thermal spectral density at fixed | q | can be 
expressed as: 

p{q,T) = i /mnf-*(g,r) tanh (^^^ . (8) 

To proceed, we need to know the fermion propagator at finite temperature. The thermal fermion propagator at real 
time is given as: 

Sik) = (7,, fc^ + mg) [ -2nn (|fco|) <5(fc' - ) ) , (9) 



here, n(x) is Fermi distribution function, 



n{x) — [exp{f3x) + 1] 



(10) 



Using the above propagator, we find the foUowing expression for the imaginary part of the correlation function at 
I q 1= hmit: 



(11) 



where, uj = xl'i^'o + ^ ■ After some straightforward calculations, the thermal spectral density is obtained as: 



1 



1 - 2 n 



(12) 



sttm ' * y V V 2 

In the non-perturbative part, the main contribution comes from the two gluon condensates since the heavy quark 







FIG. 1. (a): Bare loop diagram (b): Diagrams corresponding to gluon condensates. 

condensates are suppressed by inverse powers of the heavy quark mass and can be safely removed. The gluon 
condensate diagrams are represented in part (b) of figure (UJ. In order to calculate nonperturbative contributions, we 
use Fock-Schwinger gauge, x^Al^^ix) = 0. In momentum space, the vacuum gluon field is expressed as: 



d 



and in calculations, we use the quark-gluon-quark vertex as: 



(13) 



(14) 



where k' is the gluon momentum. 

After straightforward calculations, the non-perturbative part in momentum space is obtained as: 



Y^nonpert 

dx 



:{3(a. 



^(-1 + x)^{-l - 2a; + 2x^) - \2m^Qq^x{-l + xf 



'""2887r(TO| + g2 a; (-l + .x))4 
x(l + 4a; - 12a;^ + Qx^) + m^g2(-15 + 156a; - 441a;^ + 434a;^ - 134x*) - as(ii"e^^w'^) Aq^{-1 + x)[q^{-l + xf 
xa;2(9 + llx - Ux^ + 12a;^) + m^(-15 + 135x - 246a;^ + 176x^) + 4m^(j^a;(-3 - 8.t + 28a;^ - Mx^ + llx^)^ 
-lQ{-l + x){q.uf{q'^x^{~l + xf{'d + llx- Ux^ + 12a;^) + to^(-15 + 135x - 246a;2 + 176a;^) 
+4m^g2a;(-3 - 8a; + 28a;2 - 34a;3 + ITa;"*)) |, (15) 
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where, four- vector u'^ is the velocity of the heat bath and it is introduced to restore Lorentz invariance formally in the 
thermal field theory. In the rest frame of the heat bath, = (1, 0, 0, 0) and v? — \. In deriving the above expression, 
we have used the following relation considering the Lorentz covariance [13] : 

{Tr^GapGxa) = {gaXgpa - gaagi3\)A - {UaUxg^a - MaU<T5/3A " UpUxgaa + Ui3Uagax)B. (16) 

Contracting indices on both sides, we obtain, 

A = YA^GlpC^-P) + \{u-&%uP), (17) 

5= (18) 
where, 9^^ is the traceless, gluonic part of the stress-tensor of the QCD and it is defined as: 

= -GlxGf + \g^pGl„G^^\ (19) 

Matching the phenomenological and QCD sides of the correlation function, the sum rules for the mass and decay 
constant of scalar meson are obtained. To suppress the contribution of the higher states and continuum, Borel 
transformation over as well as continuum subtraction are performed. As a result of the above procedure, we obtain 
the following sum rule for the decay constant: 

m|(T)/|(r)e^^ = W ds p{s) e^IT^ + i?n"°"P-* , (20) 
yJi^nl J 

where A'P is the Borel mass parameter and so{T) is the temperature dependent continuum threshold. The sum rules 
for the mass is obtained applying derivative with respect to — -pr to the both sides of the sum rule for the decay 
constant of the scalar meson in Eq. pO| and dividing by itself: 

ds s p{s) exp{-^) + nr"^^^*(M2, T) 

mliT) = m ' (21) 

/4™f ds p{s) exp{^^) + Bn~-*(M2,T) 

where 

^nonpert ^^2 ^ ^ _ __i_Sn"°"P^'-* {M^ , T) , (22) 

and ijn"°"P'^''* (M^ , T) shows contribution of the gluon condensates in Borel transformed scheme and it is given by: 



+ 9M^{-1 + xfx^{-l -2x + 2x^) - 3to^M'*(-1 -f xfx^{~h + lx- Ylx^ + <6x^) 



+ 2TO^Ar j;(3 - 21a; + A&x^ - 41a;^ + \\x^) j - 4as(e9) (^TOq(1 - 2x)^(-3 + 5a;) 

+ + xfx^{9 + 11a: - 14a;^ + 12a;^) + 2TO^M2a;(3 - 23a; + SSa;^ - hlx^ + 19a;^) 



m\ M'^ (-1 + xfx^{-\h + 11a; - 26a;2 + 32a;3) ^ | 



(23) 



We use the gluonic part of energy density both obtained from lattice QCD [48| and chiral perturbation theory [4£| . 
In the rest frame of the heat bath, the results obtained in [48] at lattice QCD are reproduced well by the following fit 
parametrization for the thermal average of total energy density, (0), 

(9) = 2(69) 6 X 10-6ea;p[80(r - 0.1)](Gey^), (24) 

where temperature T is measured in units of GeV and this parametrization is valid in the interval 0.1 GeV < T < 
0.17 GeV . Note that the total energy density has been known for T > in the chira.1 perturbation theory, while this 
quantity has only been calculated for T > lOOMeT^ in lattice QCD (see [s^l and [l^l). In low temperature chiral 
perturbation limit, the results presented in [4^ are better described by the expression, 

(6) = (ej:) + 3 p, (25) 
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where, p is pressure and (OJ^) 



III! 

270 FS 



In 



is trace of the total energy momentum tensor. They are given as: 

T 



Ap 
T 



p^3T 



V 2 TT 



1 



15 T 105 

+ 



8 128 to2. 







^ exp 


T . 



(26) 



where Ap = 0.275GeV, = 0.093GeV" and = O.UGeV. 

Our final task in this section is to introduce the temperature dependent continuum threshold, so(T), gluon con- 
densate, (G^) and the strong coupling constant. The temperature dependent continuum threshold [50] and gluon 
condensate [H, are well described by the following fit parameterizations: 



s(T) = So 



4 mi 



(27) 



(G2 



IG^IO) 



12 



(4 



1.05 



1 



(28) 



where T* = 1.1 x Tc — 0.176Gey, and sq and (0|G^|0) are the continuum threshold and the gluon condensate in 
vacuum, respectively. These parameterizations are valid only in the interval < T < 170 MeV. Here, we should 
stress that the continuum threshold presented above is equal to the continuum threshold in vacuum at T = but 
it starts to diminish increasing the temperature such that at T = T* it reaches the perturbative QCD threshold, 
4mQ. This parametrization belongs to heavy- heavy system and differ considerably with the case of light-light and 
heayy-light quark systems, where the continuum threshold is related to the thermal light quark condensate (for details 
see [50l|). 

We also use temperature dependent strong coupling constant [Ell [13] as: 
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, 11 , /27rr\ 



51 



In 



MS 



, / 27rT \ 
21n(^ j 



VA- 



MS 



(29) 



where, Ajj-g ~ rc/1-14 and in numerical calculations, instead of the as in front of (0^) in Eq. (|23p the a{T) 



2.096 aP^''*(T) has been used, where aP''''\T) 



4ir 



(for details see (Slf). 



III. NUMERICAL ANALYSIS 



Present section is devoted to the numerical analysis of the sum rules for the mass and decay constant of the 
heavy scalar mesons. In further analysis, we use rric = (1.3 ± 0.05)GeT^, = (4.7 ± 0.1)Gey and (0 ] ^a^G^ ] 
0) = (0.012 ± 0.004)GeF^. The sum rules for the mass and decay constant also contain two auxiliary parameters, 
continuum threshold sq and Borel mass parameter Af^. The standard criteria in QCD sum rules is that the physical 
quantities should be independent of the auxiliary parameters. However, the continuum threshold sq is not completely 
arbitrary but is related to the energy of the first exited state with the same quantum numbers and can depend on the 
Borel mass parameter [ssl ]. Therefore, the standard criteria, does not render realistic errors, and in fact the existent 
error should be large. Hence, we will add also the systematic errors to the numerical results. We choose the values 
So = (110 ±4) GeV'^ and sq = (18 ±2) GeV'^ for the continuum threshold at Xbo and Xco channels, respectively. The 
working region for the Borel mass parameter is determined requiring that not only the higher state and continuum 
contributions are suppressed but also the contribution of the highest order operator should be small, i.e., the sum 
rules for the mass and decay constant should converge. As a result of the above procedure, the working region for 
the Borel parameter is found to be 8 GeV"^ < < 20 GeV^ for XcO and 15 GeV^ < < 30 GeV^ for Xbo mesons. 
The dependences of the masses and decay constants at T = on Borel mass parameter are shown in Figs. [211S1 These 
figures depict that the observables depend very weakly on the Borel mass parameter, Af ^ in the working regions. 

The dependence of the mass and decay constant of the Xbo and XcO mesons on temperature are presented in Figs. 
[6][9l In these figures, we show the results obtained using both lattice QCD and chiral perturbation limit values for the 
gluonic part of energy density. These figures depict that the results depend very weakly on the values of the gluonic 
part of energy density, i.e., both values obtained from lattice and chiral limit have approximately the same predictions 
in the interval, 0.1 GeV <T< 0.17 GeV at which the lattice results are valid. These figures also show that the masses 
and decay constants don't change up to T ~ 100 MeV but they start to diminish with increasing the temperature 
after this point. At near the critic or deconfinement temperature, the decay constants reach approximately to 25% of 
their values in vacuum, while the masses are decreased about 6% and 23% for bottom and charm cases, respectively. 
From these figures we deduce the results on the decay constant and mass in vacuum as presented in Tables I and II. 
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fx.oiMeV) 


Present Work 


175 ± 55 


343 ± 112 


QCD sum rules [1] 




359 


Cornell potential model [54] 




338 


QCD sum rules [55] 




510 ± 40 



TABLE I. Values of the leptonic decay constants of the heavy scalar, xto and XcO mesons in vacuum. These results have been 
obtained using the values sq = 110 GeV^ and = 17 GeV^ for Xbo, and so = 18 GeV^ and — 9 GeV^ for XcO mesons. 









Present Work 


10.10 ± 1.75 


3.71 ±0.62 


Experiment [IT] 


9.85944 ± 0.00042 


3.41475 ± 0.00031 



TABLE IL Values of the mass of the heavy scalar, Xbo and XcO mesons in vacuum. 



The quoted errors in these Tables are due to the errors in variation of the continuum threshold, Borel mass parameter 
and errors in other input parameters as well as the systematic uncertainties. Table I also include a comparison of the 
decay constant of charm case with the existing predictions of the same framework or other nonperturvative approaches. 
From this Table, we see that our predictions on the decay constant of thevco zero temperature are well consistent 
with the predictions of QCD sum rules jlj] and Cornell potential model [54] predictions, but differ considerably from 
the result obtained in [55] when the central values are considered. In Table II, we also compare our predictions on the 
masses of the heavy scalar mesons with the existing experimental data which are in a good consistency. Our results 
for the leptonic decay constants as well as their behavior with respect to the temperature can be verified in the future 
experiments. 
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FIG. 2. The dependence of the decay constant of scalar Xbo meson on the Borel parameter, in vacuum at three fixed values 
of the continuum threshold. 
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FIG. 3. The same as FIG. 2 but for mass of the scalar X6o meson. 
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FIG. 4. The dependence of the decay constant of scalar Xco nieson on the Borel parameter, in vacuum at three fixed values 
of the continuum threshold. 
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FIG. 5. The same as FIG. 4 but for mass of the scalar XcO meson. 



11 




0.00 0.04 0.08 0.12 0.16 

T(GeV) 



0.20 



FIG. 6. The dependence of the decay constant of scalar Xbo meson on temperature obtained using both lattice QCD and chiral 
perturbation limit values for the gluonic part of energy density. The values so = 110 GeV^ , and — 17 GeV^ have been 
used for the continuum threshold and Borel mass parameter in vacuum, respectively. 
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FIG. 7. The same as FIG. 6 but for mass of the scalar Xbo meson. 
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FIG. 8. The dependence of the decay constant of scalar XcO meson on temperature obtained using both lattice QCD and chiral 
perturbation limit values for the gluonic part of energy density. The values so — 18 GeV^, and = 9 GeV^ have been used 
for the continuum threshold and Borel mass parameter in vacuum, respectively. 
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FIG. 9. The same as FIG. 8 but for mass of the scalar Xco meson. 



